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SECTION  1 


INTRODUCTION 


The  properties  of  an  atmosphere  can  be  investigated  by  direct 
methods  and  indirect  methods.  The  direct  methods  have  the  advan¬ 
tage  of  providing  accurate  point  measurements,  but  the  sampling  of 
large  volumes  is  difficult  or  impossible.  The  indirect  methods, 
passive  and  active,  permit  probing  of  very  large  volumes  with  loss 
of  point  accuracy.  There  are,  of  course,  cases  where  only  one  or 
neither  of  the  methods  can  be  usefully  employed. 

The  Event  Dial  Pack  of  23  July  1970  provides  a  specific  ex¬ 
ample  of  a  situation  where  indirect  methods  of  prooing  can  be 
gainfully  employed.  During  the  early  phases  of  such  an  event,  it 
is  extremely  difficult  to  obtain  in  situ  measurements  while  both 
radar  and  microwave  radiometer  techniques  yield  useful  data, 
though  they  require  careful  interpretation.  This  report  presents 
the  results  of  a  program  directed  toward  analyzing  the  microwave 
radiometric  data  obtained  during  Event  Dial  Pack. 

Additional  micrcwave  radiometric  data  was  taken  during  this 
program  in  a  controlled  environment  to  support  the  analysis.  Both 
active  and  passive  data  were  analyzed  as  this  event  represents  the 
only  known  experiment  its  type  when  both  techniques  can  be 
compared . 


TECHNICAL  BACKGROUND 


1.  GENERAL 

The  microwave  emission  of  the  clear  atmosphere  viewed  by  a 
radiometer  from  the  ground  is  unpolarized.  When  the  atmosphere 
becomes  turbid  due  to  the  presence  of  suspended  material,  observed 
radiation  intensities  increase  and  possibly  also  exhibit  a  degree 
of  polarization.  The  cause  of  the  polarization  is  the  presence, 
in  sufficient  numbers,  of  particles  with  high  scattering  effi¬ 
ciencies,  and  observed  polarization  can  be  the  result  of  the 
scattering  process  and/or  a  property  of  the  radiation  impinging 
on  these  particles.  It  is  possible  from  spectral  measurements  of 
cross  polarized  intensity  to  determine  the  properties  of  the 
suspended  matter  when  illuminating  intensities  are  specified. 

2.  THE  MEASUREMENTS 

The  measurements  performed  during  Event  Dial  Pack  consist  of 
the  time  histories  of  the  horizontal  and  vertical  polarized  com¬ 
ponents  of  the  dust  cloud's  radiation  field  at  10.2  and  30  OH?. 
During  the  detonation  and  cloud  rise  phase  of  the  observations. 
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the  radiometers  were  directed  such  that  the  line-of-sight  passed 
directly  over  ground  zero.  Following  detonation,  the  dust  cloud 
ascended  through  the  radiometer  field  of  view.  During  the  second 
phase  the  elevation  angle  was  increased  to  48,5  degrees  to  observe 
the  overhead  passage  of  the  cloud.  A  boresighted  camera  provided 
pointing  information  for  interpretive  purposes, 

3.  THE  PROBLEM 

The  specific  problem  considered  during  the  experiment  was 
the  determination  of  physical  properties  of  a  near-earth  dust 
cloud  from  microwave  radiometric  observation.  Unfolding  the  total 
time-dependent  radiation  field  to  yield  the  properties  of  the 
cloud  was  complicated  by  the  temporal  microphysical  and  macro¬ 
physical  properties  of  the  cloud,  pertinent  geometry,  and  the 
characteristics  of  the  radiometers.  In  addition,  the  spatial 
distribution  of  natural  sources  of  radiation,  i.e.,  the  ground, 
atmosphere,  and  sun,  played  a  role  in  the  thermalization  of  the 
cloud  as  well  as  being  contributing  sources  if  scattering  is  an 
important  mechanism.  Short  of  developing  a  complete  and  general 
solution  to  this  involved  problem,  approximations  and  idealiza¬ 
tions  can  be  employed  in  models  of  the  physical  situation  to 
isolate  important  mechanisms,  provide  an  explanation  for  observed 
intensities  and  polarizations,  and  to  recover  certain  properties 
of  the  cloud.  A  goal  of  eventually  obtaining  the  full  formal 
solution  was  in  iuiiiu  during  the  program. 

4.  STATEMENT  OF  WORX 

The  proposed  tasks  for  this  program  fall  into  three  major 
categories  -  namely  microwave  data  analysis,  scattering  phenomena 
analysis,  and  passive-active  signal  comparisons.  Specific  tasks 
to  bo  performed  are  outlined  as  follows: 

a.  Microwave  data  analysis 

(1)  Reduce  radiometer  voltages  to  equivalent  trans¬ 
mission  losses. 

(2)  Estimate  beamfill  ratio  of  the  cloud  versus  time. 

(3)  Calculate  the  normalized  insertion  loss  from  cloud 
dimensions  and  positions  as  a  function  of  time. 

(4)  Conduct  a  set  of  simple  tests  to  measure  small 
particles  with  a  microwave  radiometer  to  determine 
their  signature  under  laboratory  conditions. 

(5)  Measure  the  dielectric  constant  of  actual  test 
samples . 
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b.  Scattering  Phenomena 

(1)  Determine  the  primary  mechanisms  responsible  for 
the  observed  dust  cloud  microwave  brightness  tem¬ 
peratures  and  polarizations. 

(2)  Formulate  the  mathematical  model  relating  the 
observed  brightness  temperatures  to  the  test 
geometry  in  terms  of  the  primary  radiation  transfer 
mechanism . 

(3)  Analyze  the  existing  radiometric  measurements  in 
accordance  with  the  mathematical  model  to  determine 
the  dust  cloud  parameters. 

(4)  Recommend  procedures  for  optimizing  data  return  from 
microwave  radiometric  measurements  in  future  experi¬ 
ments  . 


c.  Passive-Active  Signal  Comparisons 

(1)  Compare  radar  (X  and  Ka  band)  returns  with  radio¬ 
metric  signals. 


Calculate  back  scattering  signals  to  fcc  expected  from 
cloud  models  obtained  above.  Iterate  model  and  meas¬ 
ured  returns  to  refine  microwa /e  models 


(3)  Calculate  complete  scattering  from  dust  cloud. 


The  discussions  which  follow  do  not  necessarily  follow  the 
outlined  format  in  the  Statement  of  Work  as  the  tasks  are  cate¬ 
gory  groupings  only.  The  microwave  experimental  data  collected 
during  the  program  was  specifically  oriented  toward  supporting 
the  analysis  of  Event  Dial  Pack  data. 


bo  major  difficulties  were  encountered  in  performing  all 
tasks  except  for  differentiating  cloud  thermal  changes  from  in¬ 
sertion  loss  changes.  The  radar  data  was  very  minimal  for  a  full 
comparative  analysis. 

5.  MEASUREMENTS 

a.  General 


The  measurements  reported  upon  herein  were  conducted  to 
provide  ancillary  data  to  that  obtained  during  Event 
Dial  Pack.  All  experiments  were  conducted  under  con¬ 
trolled  conditions  for  the  purpose  of  determining  atten- 
uauxunf  I ai i z a 1 1 o n  and  refractive  properties  on  duet 

particles.  These  measurements  were  necessary  to  assure 
that  the  analysis  was  approached  in  a  meaningful  way. 
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The  data  was  further  reduced  during  the  program  to  pro¬ 
vide  beamfill  data  for  the  microwave  radiometer  artennas. 
These  data  along  with  dielectric  constant  data  are  re¬ 
quired  for  an  evaluation  of  insertion  loss  and  trans¬ 
mission  loss. 

Experimental  Data 

(1)  Transmission  Loss  Data 

Initially  the  experiments  of  this  program  were 
planned  as  scaled  simulations  of  Event  Dial  Pack 
dust  clouds.  The  clouds  were  to  be  created  at  close 
range  with  a  constant  density  and  dust  particle  size. 
Subsequent  post  contract  evaluation  indicated  that 
the  proposed  approach  was  too  difficult  to  instru¬ 
ment  with  any  accuracy  and  that  cloud  density  could 
not  be  controlled  for  100%  antenna  beamfill. 

An  alternate  solution  was  used  which  proved  to  be 
satisfactory,  in  which  density  could  be  simulated 
and  easily  controlled.  Instead  of  creating  dust 
clouds,  sheets  of  dry  foam  were  placed  directly  in 
front  of  the  antennas  with  the  antennas  aimed  at  the 
zenith.  Sand  was  weighud,  measured  and  evenly 
spread  over  the  foam  in  the  antennas'  field  cf  view. 

Some  uncertainties  exist  with  this  technique,  in  that 
the  dust  particles  are  now  in  the  near  field  of  the 
antenna.  Experiments  with  larger  targets  show  that 
the  effects,  if  any,  on  the  data  are  small  as  the 
signals  measured  are  noncoherent  noise.  Uncontrolled 
differencies  in  the  simulation  are  sidelobe  and 
backlobe  contributions. 

Curves  of  differential  temperature  versus  ounces  of 
sand  are  shown  in  Figures  1  and  2.  These  data  rep¬ 
resent  both  vertical  and  horizontal  polarized  signals 
as  no  differences  were  noted  after  corrections  were 
made  for  gain  differences.  The  sinusoidal  nature  of 
the  curves  should  be  noted.  This  characteristic  is 
typical  of  abrupt  interfaces  between  materials  of 
differing  dielectric  constants.  Each  complete  cycle 
represents  one-half  wavelength  of  material  thickness. 
Damping  of  the  sinusoid  represents  the>  losses  in  the 
material.  For  large  dust  clouds,  this  effect  should 
disappear. 

(2)  Dielectric  Measurements 

The  dielectric  constants  of  the  materials  used  during 
this  program  were  measured  using  a  large  coaxial  line 
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10.2  GHz  Differential  Radiometric  Temperature 


terminated  in  50  ohms.  Materials  were  placed  in  the 
dielectric  region  of  the  coaxial  line  and  reflection 
coefficients  were  measured  using  a  Time  Domain  Re- 
flectometer  (TDR)  technique.  The  dielectric  con¬ 
stants  for  the  fine  and  coarse  sands  used  were  2,64 
and  3.51,  respectively.  Samples  of  fine  and  coarse 
Canadian  soil  were  obtained  after  the  measurements 
portion  of  this  program;  however,  dielectric  measure¬ 
ments  were  made  in  the  same  manner  and  found  to  be 
3.86  and  4.40  for  the  fine  and  coarse  soil,  respec¬ 
tively. 

(3)  Antenna  Beamfill  Ratio 

The  only  means  for  evaluating  antenna  oeamfill  ratio 
during  Event  Dial  Pack  was  from  the  16mm  movie  bore- 
sight  camera  film.  The  camera  had  a  25mm  lens  and 
microwave  radiometer  antenna  half-power  (3db)  beam- 
widths  were  4.0  degrees  at  10.2  GHz  and  2.4  degrees 
at  30  GHz.  Thus,  the  antenna  beam's  intercept  area 
could  be  calculated  assuming  the  dust  cloud  distance 
to  be  equal  to  ground  zero.  An  overlay  was  con¬ 
structed  equal  to  this  area  and  percentage  beamfill 
was  measured  from  the  film.  A  plot  of  beamfill  ratio 
calculated  by  this  technique  is  shown  in  Figure  3. 
Beamfill  ratios  after  panning  to  the  cloud  were 
always  100%. 
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SECTION  II 


ANALYSIS 


1.  GENERAL 

The  approach  taken  during  the  analysis  was  to  first  examine 
the  radiative  properties  of  dust  particle  distributions  obtained 
by  aircraft  sampling  at  times  exceeding  10  minutes  after  detona¬ 
tion.  As  earlier  aircraft  data  was  not  available,  a  variety  of 
extrapolations  to  earlier  times  were  made  in  an  attempt  to  simu¬ 
late  early  phase  radiometric  observations. 

A  second  approach,  whereby  particle  refractive  index  and  an 
arbitrary  particle  size  distribution  would  have  been  subjected  to 
iterative  adjustment  to  match  emission,  polarization  and  radar 
backscattering  observations,  was  abandoned  for  lack  of  reliable 
data.  This  second  approach  would  have  provided  an  assessment  of 
the  validity  of  the  extrapolation  of  physical  sampling  data  along 
with  an  alternate  determination  of  cloud  radiation  properties. 

2.  QUALITATIVE  DATA  INTERPRETATION 

The  data  acquired  during  the  Event  Dial  Pack  is  shown  in 
Figure  4.  Time  progresses  to  the  left  and  all  significant  events 
are  labeled.  The  outputs  of  the  horizontally  and  vertically 
polarized  receiver  channels,  designated  H  and  V,  respectively, 
are  displayed  for  both  10.2  GHz  and  30  GHz.  The  differences,  H-V, 
which  are  proportional  to  polarization  are  also  displayed.  Cali¬ 
brations  for  all  channels  are  1°C  per  millivolt,  giving  25°C  full 
scale  for  the  H  and  V  channels  and  5°C  full  scale  for  the  polar¬ 
ization  difference  channels. 

The  results  for  10.2  GHz  radiometer  are  shown  in  the  top 
three  traces.  These  data  are  contaminated  by  signals  from  a 
nearby  radar,  which  prior  to  the  event  were  directed  away  from 
ground  zero  and  toward  the  radiometer  location.  At  TQ,  this  radar 
was  pointed  toward  ground  zero,  but  some  interference  is  still 
evident.  As  the  level  of  this  interference  does  not  appear  to  be 
time  dependent,  it  is  concluded  that  radar  energy  scattered  from 
the  cloud  toward  the  radiometers  is  not  detectable. 

The  behavior  of  the  polarization  channel,  although  irregular 
and  indicating  some  degree  of  polarization,  was  difficult  to  in¬ 
terpret.  Figure  5  gives  boresight  camera  frames  corresponding  to 
the  times  labeled  1-4  on  the  strip  chart  record.  Thus,  at  time  2, 
the  radiometers  were  boresighted  on  the  cloud,  and  as  the  geometry 
of  Figure  6  depicts,  the  radiometer  beams  were  completely  filled. 
The  H  and  V  channels  follow  each  other  quite  closely  from  the  be¬ 
ginning  of  the  cloud  rise  through  time  4.  A  notable  character¬ 
istic  of  these  records  is  that  the  polarization  channel  exhibited 


9 


Figure 


-  10  GHz  and  30  GHz  Radiometric  Response1 
(View  Angle  =  17°  above  the  horizon) 
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Figure  6  -  Diagram  of  Experiment  Geometry 


its  largest  excursions  when  the  rate  of  changa  of  H  and  V  were 
greatest.  The  conclusion  is  that  small  differences  in  gain,  time 
constant,  H  and  V  plane  beam  shapes,  and  nonlinearity  are  the 
most  probable  causes  of  the  apparent  small  and  variable  degree  of 
polarization. 

The  results  from  the  30  GHz  radiometer  are  shown  in  the 
bottom  three  traces  of  Figure  4.  If  the  radiation  received  from 
the  direction  of  the  cloud  truly  possess  a  degree  of  polarization, 
it  should  be  evident  in  this  data. 

The  immediate  striking  feature  of  the  30  GHz  records  is  the 
large  change  in  the  polarization  occurring  when  the  rising  cloud 
passes  through  the  radiometer  field  of  view.  Otherwise,  the 
polarization  channel  exhibits  little  variability  and  substantiates 
the  conclusion  arrived  at  from  examination  of  the  10.2  GHz  data. 

Closer  scrutiny  of  the  H  and  V  traces  shows  that  while  they 
follow  each  other  closely,  the  peak-to-peak  excursions  of  the 
larger  scale  structure  shows  evidence  of  differential  gain.  As 
there  appears  to  be  no  radar  interference,  the  ratio  of  the 
peak-to-peak  excursions  on  the  two  channels  can  be  scaled  for  the 
short  period  just  prior  to  TQ.  During  this  period,  clear  sky  and 
broken  cumulus  clouds  were  in  the  field  of  view  so  that  polariza¬ 
tion  would  not  be  expected.  The  ratio  AH/AV  is  approximately 
0.67.  The  ratio  of  the  changes  in  H  and  v  due  to  the  dust  cloud 
at  time  2  is  also  approximately  0.67.  From  this,  it  was  concluded 
that  differential  gain  existed  and  that  the  vertical  channel  had 
the  highest  gain.  If  the  polarization  channel  is  corrected  for 
the  indicated  amount  of  differential  gain,  the  amplitude  of  its 
variations  over  the  entire  record  is  reduced  to  noise. 

The  true  degree  of  polarization,  ir,  would  be  less  than  or 
equal  to  (H-V)/(H+V).  Therefore,  if  the  absolute  base  radiation 
levels  were  assumed  to  be  zero,  and  if  differential  gain  were 
ignored,  the  degree  of  polarization  at  time  2  would  be 

"  £  SSTTTF  =  °-20 

or  approximately  20  percent.  To  be  more  realistic,  however,  a 
background  level  of  40°K,  which  corresponds  to  a  base  level  of 
400  mv,  can  be  assumed  (see  Figure  7).  Under  this  assumption 

45 

n  ^  4£6'+33fe  =  0,0439 

or  approximately  4.4  percent.  This  would  be  a  measurable  and  sig¬ 
nificant  amount  of  polarization.  However,  when  differential  gain 
is  accounted  for,  the  degree  of  polarization  becomes  inconsequen¬ 
tial. 
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Figure  7  -  Sky  Brightness  Temperatures  for 
A  Clear  Standard  Atmosphere 


1  4 


Despite  the  conclusions  of  this  qualitative  look  at  the  data, 
the  problem  was  further  examined,  using  all  available  data,  first 
with  Rayleigh  calculations  and  then  with  extensions  to  Mie  theory. 
After  this,  a  final  assessment  of  the  microwave  radiometric  data 
and  its  significance  was  made. 

3.  RAYLEIGH  APPROXIMATIONS  -  LATF  PHASE 


As  a  starting  point  to  quantitative  analysis,  the  aircraft 
particle  sampling  data  for  Event  Dial  Pack,  pass  #1,*  T_  4  11 
minutes,  are  examined  and  then  used  as  the  basis  for  calculations. 
The  particle  size  distribution  is  given  in  Table  I,  For  this 
distribution  of  669  particles,  the  average  particle  size  (dia¬ 
meter)  is  59  microns r  giving  a  circumference  to  wavelength  ratio 
(x  =  ird/X)  of  0.018  at  30  GHz,  well  within  the  Rayleigh  criteria. 
Even  the  largest  particles  in  the  distribution  (300y)  have  the 
size  parameter  x  <  0.1.  It  is  therefore  assumed  that  the  Rayleigh 
approximations  for  extinction  and  absorption  are  valid. 


a.  The  Equivalent  Volume  Radius 

If  all  of  the  particles  in  a  size  distribution  are  with¬ 
in  the  Rayleigh  range,  then  the  true  distribution  can  be 
replaced  for  computational  purposes  by  the  same  total 
number  of  particles  of  equivalent  volume  radius,  re.  By 
this  procedure,  the  extincu^an  cross-sccticn,  and 

the  volume  extinction  coefficient,  6ex,  can  be  simply 
calculated  by  algebraic  means. 


The  volume,  VD,  of  the  669  dust  particle  is  given  by  the 
summation 


I 


D(ri)r? 


(1) 


where  Dtr^)  is  the  tabulated  dust  particle  size  distri¬ 
bution  of  Table  I  and  ri  is  particle  radius.  The  sub¬ 
script  i  has  been  introduced  to  indicate  the  discrete 
sample  values.  Performing  this  numerical  integration 
produces  the  following  table  of  values  (Table  II)  with 
the  result  that  Vp  =  3.75x10®  y3.  Equating  V'd  to  the 
volume  of  669  particles  of  radius  re  and  then  solving 
for  re  gives 

re  =  5.1  y 


^Report  uAsA  z600,  Mk1/0,  *'K-y4B 
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TABLE  I 


DUST 

PARTICLE  SIZE  DISTRIBUTION  FOR 

PASS  #1* 

Diameter 

Number  of 

Diameter 

Number  of 

U 

Particles 

y 

Particles 

10 

4 

130 

17 

20 

9 

140 

9 

30 

65 

150 

14 

40 

86 

160 

10 

50 

95 

170 

9 

60 

85 

180 

4 

70 

68 

190 

7 

80 

49 

200 

9 

90 

39 

220 

1 

100 

37 

250 

5 

110 

25 

300 

1 

120 

20 

350 

1 

Mean  Par 

tide  Diameter 

5911 

Total  Number  of  Particles 

-  669 

Total  Dust  Collected  - 

13.9  ±  l.Omg 

Sampling 

Times:  On  T_  + 

10  minutes,  51 

seconds 

Off  Tq  + 

11  minutes,  32  seconds 

Total  Volume  Sampled  -  1.82ms 


Cloud  Width  -  1920m 


*  Data  taken  from  Report  DASA  2600  MR170  FR-948 
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TABLE  II 


Table  of  Values  for  Calculating  the 
Equivalent  Volume  of  the 
Pass  #1  Particle  Distribution 


i 

D(ri) 

r±  (V) 

r?(ys) 

^VD(y3) 

1 

4 

5 

1.25xl02 

2 . 09x10  3 

2 

9 

10 

1 . 00x10  3 

3.77x10** 

3 

65 

15 

3.38x10  3 

9.17x10s 

4 

86 

20 

8 . OGxlO  3 

2.88x10s 

5 

rv  r 
-/ 

25 

1.56x10 

6 . 21x10 6 

6 

85 

30 

2.70x10** 

9.61x10s 

7 

68 

35 

4.29x10** 

1 . 22x10  7 

8 

49 

40 

6.40x10** 

1 . 31x10  7 

9 

39 

45 

9. 11x10  4 

1.49x10 

10 

37 

50 

1. 25x10 ! 

1 » 94x10  7 

11 

25 

55 

1 . 66x10  J 

1.7 4x10  7 

12 

20 

60 

2 . 16x10  5 

1 . 81x10  7 

13 

17 

e  ^ 

2.75xl05 

1. 96x10  7 

14 

9 

70 

3.43x10s 

1 . 29x10  ~ 

15 

14 

75 

4.22x10s 

2.48x10 

16 

10 

80 

5.12x10s 

2 . 14x10  7 

17 

9 

65 

6.14xl0S 

2.31x10 

18 

4 

90 

7.29x10 

1.22x10 

19 

7 

95 

8,57xlOS 

2 . 51x10  7 

20 

9 

100 

1.00x10 

3 . 77x10 7 

21 

1 

110 

1.33x10s 

5.57x10 

22 

5 

125 

1.95x10 

4 , 08x10  7 

23 

1 

150 

3.38x10s 

1.4 2x10 7 

24 

1 

175 

5.36x10 

2.24x10 

VD  =  3.745x10®  y® 


This  equivalent  particle  radius  is  1.73  times  larger 
than  the  average  particle  radius  of  29. 5 y  for  the  dis¬ 
tribution  and  the  corresponding  value  of  wavelength 
ratio  x  at  30  GHz  is  0.032,  well  within  the  Rayleigh 
criteria.  Note  that  no  significance  should  be  given 
to  re  other  than  its  value  in  simplifying  calculations. 

b.  Simulation  Experiment  Results 

Changes  in  measured  apparent  temperature  were  recorded 
as  a  function  of  the  total  weight  in  ounces  of  the  sand, 
which  was  spread  uniformly  over  a  polyurethane  sheet 
(Figures  1  and  2).  This  data  as  a  function  of  mass  per 
unit  column  for  the  30  GHz  radiometer  vertical  polariza¬ 
tion  channel  for  both  grades  of  sand  is  shown  in  Figure 
8.  The  periodicities  can  again  be  noted  by  inspection 
of  individual  points.  Analyses  of  the  data  show  that 
9.8  ounces  of  coarse  sand  and  38  ounces  of  fine  sand 
were  required  to  cause  an  increase  of  13.5°C  in  apparent 
temperature,  as  observed  during  Event  Dial  Pack.  These 
weights  of  sand  correspond  to  layer  thickness  of  0.225cm 
and  1.059  cm,  respectively.  The  mass  in  a  centimeter 
square  column  was  0.387  gins  for  coarse  sand  and  1.52  guts 
for  fine  sand.  Results  for  horizontal  polarization  are 
not  given,  as  they  were  consistent  with  those  of  the 
vertical  channel,  as  was  expected. 

The  values  cited  above  may  not  be  exact  as  the  layers  of 
sand  were  in  the  near  field,  but  are  sufficient  for  the 
comparisons  to  be  made.  It  is  expected  that  sand  par¬ 
ticles  in  suspension  may  not  yield  the  same  results  and 
that  considerably  more  mass  per  unit  column  would  be 
required  for  comparable  extinction  properties.  The  ex¬ 
cess  required  mass  depends  on  the  particle  dielectric 
properties  and  size  distribution. 

At  10*2  GHz  the  experimental  data  was  more  difficult  to 
interpret  because  of  sinusoidal  interference  effects, 
but  it  appears  that  approximately  10  ounces  of  each  grade 
of  sand  gives  rise  to  the  5.5°C  increase  in  apparent 
temperature  as  observed  during  Event  Dial  Pack. 

The  modulus  of  the  dielectric  constant,  for  the 

coarse  sand  was  determined  to  be  3.51,  while  that  of  the 
fine  sand  was  2.64.  As  the  loss  tangent  can  only  be 
implied  by  the  slope  of  the  curve,  computer  model  calcu¬ 
lations  for  dust  with  the  refractive  properties  of  these 
coarse  and  fine  sands  were  made  for  imaginary  parts  of 
the  refractive  index  of  0.01,  0.005  and  0.001.  Hand 
calculations  were  made  using  approximate  formulas  for 
0.005,  i.e.,  m  *=  (e^)  -  i0.005. 
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c.  The  Volume  Extinction  Coefficient 

Using  the  equivalent  volume  radius,  re,  for  the  dust  par¬ 
ticle  distribution  of  pass  #1,  the  volume  extinction  co¬ 
efficient  for  the  dust  cloud  at  T  +  11  minutes  can  be 
calculated.  First,  the  particle  extinction  cross  section, 
°ex»  8  calculated  from1 

3ex  =  48if?'r|  (2) 


where 


v  =  real  part  of  dust  particle  refractive  index 
k  =  imaginary  part  of  dust  particle  refractive  index 
X  =  wavelength  cf  radiation 


and 


2^  -  (v2+ic2)2  +  4  ( v 2 “k 2 )  +  4 

Then,  substituting  r3  =  1.34x\0“7  cm3,  the  particle  ex¬ 
tinction  cross-sections  for  x  =  1cm  and  x  =  2.94cm,  re¬ 
spectively,  (designated  by  superscripts  k,  x)  are  found 
to  be 


m  e  V005  m  =  1.63=i0.005 

°ex(cm2)  =  2 . 37x10” 8 

a£x(cm2)  =  6 .  t  8 . 06x10” 9 

The  volume  extinction  coefficients,  Sex,  can  now  be  cal¬ 
culated  using  particle  densities.  From  pass  #1  data, 
there  were  669  particles  collected  from  1.82m3  of  cloud. 
Therefore,  a  particle  density  of  3.68xl0“4  per  cm3  is 


arrived 

are 

at  and 

the  corresponding 

extinction  coefficients 

m  =  1.88-i0.005 

m  =  1 . 63-iO . 005 

Sex 

(cm-1 ) 

=  7 . 13x10“ 1 2 

8 . 72x10“ 1 2 

Six 

(cm'1  ) 

=  2. 4 3x10” 12 

2 . 97x10” 1 2 

'penndorf,  R.  B. ,  "Scattering  and  Extinction  Coefficients  for 
Small  Absorbing  and  Non-Absorbing  Aerosols,"  J.  Opt.  Soc.  Amer., 
'voi .  52,  No.  5,  August  1962,  pp.  836—904. 
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d.  Optical  Thickness  and  Transmittance 

The  optical  thickness  of  the  cloud,  T,  is  defined  by 


3  gv  ( s )  ds 


and  if  it  is  assumed  that  6ex  is  independent  of  position, 
s,  then 

T  =  3exW  (4) 

where  W  is  the  width  of  the  cloud.  Working  under  this 
assumption  (particle  distributions  as  a  function  of  posi¬ 
tion  are  not  available),  the  optical  thickness  for  the 
respective  wavelengths,  using  the  cloud  width  of  1.92  km 
for  pass  #1  are 


m  =  1.88-i0.005  m  =  1.63-10.005 


-k  _ 


1.37x10" 


1.67x10" 


xx  = 


4.67x10" 


5.70x10' 


Transmittance  is  related  to  optical  thickness  by 
t  =  exp  {-T ) 
from  which  we  obtain 


m  =  1.88-i0.005  m  -  1.63-i0.005 


-k  _ 


1-0.00000137 


1-0.00000x67 


1—0.000000466  x— 0.00000  /u 


These  results  indicate  that  at  T  -  11  minutes,  the  cloud 
is  for  all  practical  purposes  transparent  at  wavelengths 
of  1  and  2.94  cm.  Reference  to  Figure  9  establishes  this 
to  be  the  case,  as  the  radiometers  are  obviously  seeing 
background  radiation.  There  was  no  need  to  consider  the 
relative  importance  of  scattering  and  absorption  in  this 
case,  but  since  a  comparison  of  scattering  coefficients 
is  made  at  a  later  point,  the  volume  scattering  coeffi¬ 
cient  8gc  is  next  calculated. 
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T* 

Figure  9-10  GHz  and  30  GHz  Radiometric  Response 
(View  Angle  =  48.5°  above  horizon) 


2589 
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e.  The  Equivalent  Volume-Squared  Radius 

Previously,  the  equivalent  volume  radius  was  determined 
for  use  in  calculating  the  extinction  coefficients.  How¬ 
ever,  when  using  re  in  the  Rayleigh  formula  for  scatter¬ 
ing  cross  section,  a  significant  underestimate  results. 
Experience  has  shown  that  an  alternate  radius,  the 
equivalent  volume-squared  radius  rT/v,  defined  by 


produces  accurate  values  for  the  scattering  coefficients. 

With  the  aid  of  the  following  table  of  values  (Table  III), 
the  volume-squared  for  the  pass  #1  particle  size  distri¬ 
bution,  Vq2 ,  is  found  to  be  1.88xlCT 5  p6.  Normalizing 
Vjj2  to  the  total  number  of  particles  and  solving  for  rvv 
gives 

rw  =  74  w 

This  equivalent  particle  radius  is  2.5  times  larger  than 
the  average  particle  radius  and  1.45  times  larger  than 
re,  but  is  also  well  within  the  Rayleigh  range  with  the 
corresponding  value  of  x  of  0.046.  As  with  re,  r^ 
should  be  given  no  special  significance. 

f.  The  Volume  Scattering  Coefficient 

Under  the  assumptions  that  the  cloud  is  optically  thin 
(t<<1)  and  that  single  scattering  governed  by  Rayleigh 
laws  holds,  extinction,  absorption  and  scattering  are 
related  by 

^ex  =  ^ab  +  ®sc 

The  alternative  therefore  exists  for  calculating  3§c 
directly  or  indirectly  through  the  calculation  of  Bab. 

The  first  alternative  was  chosen  here  because  of  its 
simplicity.  The  particle  scattering  cross-section,  6SC, 
can  be  calculated  using 
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TABLE  III 


Table  of  Values  for  Calculating  the 
Equivalent  Volume-Squared  of  the 
Pass  #1  Particle  Distribution 


i  D(ri)  VD2 (y3) 


D(ri) AVD2  (y6 ) 


4 

5. 23xl02 

1.10x10s 

9 

4 . 19x10  3 

1 . 58xl08 

65 

1.42x10“ 

1. 31X101 0 

86 

3. 35x10“ 

9 . 6  3x10 1 0 

95 

6.54x10“ 

4 . 07x10 1  1 

85 

1.13x10s 

1. 09xl0x  2 

68 

1.80x10s 

2 . 20x10 1 2 

49 

2.68x10s 

3 . 52x10 1 2 

39 

3.82x10s 

5 . 69x10 1 2 

37 

5.24x10s 

1 . 02x10 1 3 

25 

6.96x10s 

1 . 21x10 1 3 

2<J 

9.05x10s 

1, C4xi0 1 3 

17 

1.15xl06 

2. 24x10  1  3 

9 

1.44xl06 

1 . 86x10 1 3 

14 

1 . 77x10s 

4 . 38x10 1 3 

10 

2.15x10° 

4 . 62x10  1  3 

9 

2.57x10s 

5 . 94x10  1  3 

4 

3.06x10s 

3.74x1c1  3 

7 

3.59x10s 

9. 03x10*  3 

9 

4.19x10s 

1.58x10*“ 

1 

5.57x10s 

3.10x10* “ 

5 

8.17x10s 

3.34x10* “ 

1 

1 . 42x10  7 

2.02x10*“ 

1 

2 . 25xl07 

5.06x10* “ 

V&2  ^  18 . 80x10 1 4  y3 
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(from  Deirmendj ian 2) 


asc  -  24’ 


•2  , 

I m  - 
|m  *+2 


r4lT  s 

I"T  rwl 


(8) 


where  m  =  v-iic  »  particle  complex  index  of  refraction. 


Note  that  the  equivalent  volume-squared  radius  just 


derived  is  used  in  this  expression.  The  resulting 
scattering  cross-sections  and  volume  scattering  co 
efficients  for  the  two  measurement  wavelengths  are 

m  =  1. 68-10 . 005 

m_ 

=  1.63-10.005 

a^c (cm2) 

=  4.4  3x10“ 1 0 

2 . 68x10' 1 0 

asc  (c™2) 

=  5 . 93x10“  1  2 

3 . 59x10“ 1 2 

8^c(cm“1) 

=  1.63x10“* * 

9 . 86x10“  1  4 

Qsc  <cn'~1) 

=  2 . 18x10“  1  5 

1 . 32x10“ 1 5 

The  important  thing  to  note  is  that  the  volume  scatter¬ 
ing  coefficients  are  more  than  an  order  of  magnitude 
smaller  than  the  volume  extinction  coefficients.  For 
all  practical  purposes,  then,  extinction  is  equivalent 
to  absorption  for  this  particular  case. 


4.  CLOUD  LINE-OF-SIGHT  TRANSMITTANCE 


The  transmittance  (transmission  loss)  of  the  cloud  along  the 
axis  of  the  radiometer  beams  can  be  approximated  from  the  measured 
changes  in  brightness  temperature,  if  it  is  assumed  that  no  differ¬ 
ences  exist  in  the  magnitude  of  the  background  brightness  and 
cloud  thermodynamic  temperauura.  ±ne  first  assumption  was  neces¬ 
sary  because  absolute  radiation  temperatures  were  not  obtained. 

The  second  assumption  is  a  simplifying  assumption  introduced  for 
lack  of  information  on  the  thermalization  of  the  cloud. 

At  the  time  when  the  radiometers  were  first  boresighted  on 
the  center  of  the  cloud  (T0  +  48  seconds)  with  their  beams  com¬ 
pletely  filled,  the  change  in  apparent  temperature  or  transmission 
loss,  ATb,  over  the  atmospheric  background  can  be  expressed  as 

=  Tbbt  +  Tbc  -  Tbb  (9) 


Deirmendj ian,  D. ,  "Complete  Microwave 
Propel Lias  uf  rulyuispersed  Cloud  ana 
Corporation,  Report  R-422-PR,  December 


Scattering  and  Extinction 
xain  elements,"  Tiie  Rand 
1963. 
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where 


Tbb  =  apparent  temperature  of  the  background 

t  =  cloud  transmittance 

Tbc  =  ®PParent  temperature  of  the  cloud 

For  computational  purposes  it  was  assumed  that  the  apparent 
temperatures  of  the  background  at  10.2  and  30  GHz  are  14  and  40°K, 
respectively  (see  Figure  7) ,  and  that  the  cloud  was  thermalized 
at  290° K . 

A  third  assumption  that  scattered-in  radiation  is  negligible 
compared  to  cloud  self -emission  was  made.  This  assumption  is  a 
good  approximation  for  optically  thin  clouds  such  as  the  one  being 
considered,  thus 

Tbc  =  (1‘t}  Tc  (10) 

where  Tc  is  the  thermodynamic  temperature  of  the  cloud.  Equations 
(9)  and  (10)  can  be  combined  to  give  the  following  expression  for 
transmittance : 


ATb 

t  =  1  -  -  (11) 

Tc  “  Tbb 

Substituting  appropriate  values  in  equation  (11),  the  approximate 
cloud  transmittances  at  each  frequency  are 

tk  =  0.946 

tx  =  0.979 

To  go  one  step  further,  substitution  of  these  transmittances 
back  into  equation  (10)  yields  the  apparent  temperatures  of  the 
cloud  alone,  from  which  insertion  loss  can  be  implied 

t£c  -  15*7°K 

Tgc  =  6 • 05  °K 

The  ratios  of  the  cloud  apparent  temperatures  is  2.6,  which 
is  somewhat  less  than  the  theoretical  value  of  2.94  (the  ratio  of 
the  wavelengths)  expected  for  Rayleigh  particles,  but  with  the 
assumptions  employed  no  significance  can  be  placed  on  this  out¬ 
come  . 


Th<=>ne  results  will  be  useful  for  comparative  purposes  at  a 
later  point  when  similar  results  are  derived  from  t  e  computer 
models . 
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5.  RAYLEIGH  APPROXIMATIONS  -  EARLY  PHASE 


Any  calculations  of  cloud  properties  at  the  early  stages  of 
its  evolution  will  be  truly  speculative  as  supporting  in-situ 
measurements  are  not  available.  Interpretation  of  radar  data  in 
terms  of  particle  sizes  and  densities  should  yield  some  reason¬ 
able  results.  However,  it  is  interesting  to  derive  radiometric 
results  independently  from  extrapolations  of  the  aircraft  par¬ 
ticle  sampling  data. 


a.  The  Dust  Cloud  at  T0  +  48  Seconds 

If  the  mass  density  from  the  period  of  30  to  60  minutes 
after  blast  is  linearly  extrapolated  back  to  the  time 
T0  +  48  seconds,  it  is  found  that  the  mass  density  would 
have  been  approximately  30xlQ-8  gms  cm-3  (see  Figure  10). 
This  is  an  increase  of  mass  density  by  a  factor  of  44 
over  that  at  T0  +  11  minutes.  While  it  is  not  reason¬ 
able  to  expect  that  the  particle  number  density  has  in¬ 
creased  by  the  same  factor,  as  much  of  the  mass  at 
early  times  is  made  up  cf  larger  particles,  it  is  not 
unreasonable  to  assume  this  to  be  the  case  just  to  esti¬ 
mate  the  effects  of  a  mass  increase  where  most  particles 
remain  in  the  Rayleigh  size  range.  When  this  is  done, 
the  extinction  coefficients  as  previously  udlculatea  for 
T0  +  11  minutes  are  just  increased  by  a  factor  of  44. 

That  is, 

m  =  1. 88-iO. 005  m  =  1. 63-iO. 005 
aL^nr1)  =  3.14xlO-1  0  3 . 84x10“ 1 0 

vA 

=  1 . 07x10“ 1  0  1.31x10-'° 

From  Figure  11  the  cloud  width  at  T0  +  48  seconds  is 
found  to  be  1.0  km  for  which  the  optical  thicknesses  and 
corresponding  transmittances  are  found  to  be 


m  =  1. 88-iO . 005 
3.14X10'5 
1.07X10-5 
1 
1 


m  =  1 . 6  3-iQ .0  05 
3 . 84x10“ J 
1.31xl0“5 
1 
1 


These  results  indicate  that  the  linear  number  density 

.  i —  . n  i  _  r.  j -  n 3  *-  <■  ■ 

uruii  axuac  win  ^iuuubc:  uic  mcasuiau 

and  that  a  much  larger  density  would  be  required  to  give 
cloud  apparent  temperatures  in  the  range  of  those  od- 
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Figure  10  -  Linear  Extrapolation  of  Mass 
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Figure  11  -  Early  Phase  Cloud  Dimensions 


{Data  taken  from  DASA  2600  MR170  FR-948) 


served  during  Event  Dial  Pack.  From  the  transmittance 
values  it  can  readily  be  determined  that  a  mass  density 
over  1000  times  higher  would  be  necessary. 

An  additional  increase  in  mass  density  by  a  factor  of 
1000  can  be  simply  achieved  by  increasing  the  equivalent 
volume  radius  by  a  factor  of  10.  This  would  be  re=510p, 
still  in  the  Rayleigh  range,  but  of  course  an  associated 
particle  size  distribution  would  have  some  particles  in 
the  Mie  range.  As  enters  into  the  expression  for  par¬ 
ticle  extinction  cross  section,  equation  (2) ,  the  factor 
of  1000  is  realized.  The  values  of  extinction  cross  sec¬ 
tion,  volume  extinction  coefficient,  optical  thickness 
and  transmittance  calculated  for  this  new  distribution 
are 


m  =  1.88-i0.005  m  =  1.63-i0.005 


Oex (cm2) 

C 

1.94x10" 5 

2. 37x10" s 

Oex (cm2) 

=5 

6 . 60x10“ 6 

8 . 06x10“ 6 

Bex (cm-1 ) 

3 . 14x10" 7 

3 . 84x10" 7 

Sex (cm-  1 ) 

= 

1.07x10" 7 

1. 31x10" 7 

Tk 

- 

3 . 14x10" 2 

3 . 84x10" 2 

TX 

D 

1 . 07x10" 2 

1.31x10" 2 

ft 

= 

0.969 

0.962 

tX 

a 

0.989 

0.987 

b.  Absorption  vs.  Scattering 

The  relative  values  of  absorption  and  scattering  are 
dependent  upon  the  dielectric  properties  of  the  particles 
as  well  as  their  size.  If  the  absorption  coefficients 
are  an  order  of  magnitude  or  greater  than  the  scattering 
coefficients,  then  the  emission  of  the  cloud  can  be 
simply  and  directly  calculated.  On  the  other  hand,  if 
the  scattering  coefficients  are  of  the  same  magnitude  or 
greater  than  the  absorption  coefficients,  then  an  itera¬ 
tive  solution  of  the  radiative  transfer  equations  is 
called  for.  For  this  reason,  it  is  essential  to  compare 
these  coefficients  for  each  particle  distribution  being 
considered. 

Using  equation  (8)  with  a  value  of  737y  for  rw ,  ten 

1.  1  —  A  «  1.U  £  a.  mm  JLU  a  mm  tm  **  il  1  V\  (  I  4-  1  M  4-  Vs  4-  />  1  1  1  V»  r-4 

UXIUCO  U44U  W  iV4.  U1»0  t'MWW  fl  *  >-•**  *-*^«  «* 

volume  scattering  coefficients  are  obtained: 


30 


1 .63-iO . 005 


3sc  (cm"  X) 


gsc (cn“  1 


)  = 


m  -  1, 88-iO.  005  rr  = 

7.13x10“ 9  4 .31x10“ 9 

9 . 55x10“  1  °  5. 77x10"  1 0 


From  these  values,  the  volume  absorption  coefficients 
are  obtained  by  taking  the  difference  of  the  volume  ex¬ 
tinction  and  scattering  coefficients,  in  accordance  with 
equation  (7).  The  results  arei 


m  =  1.88-i0,005  m  =  1. 63-iO. 005 


6ab(cm_1)  “ 

(cm"1)  = 


2.43xl0‘ 


1.06x10' 


3.41x10" 7 
1.30x10” 7 


The  effects  of  wavelength,  real  part  of  the  refractive 
index,  and  particle  size  are  now  apparent.  At  1  cm  the 
scattering  coefficient  is  within  an  order  of  magnitude 
cf  the  absorption  coefficient  for  both  values  of  refrac¬ 
tive  index,  while  at  3  cm,  it  is  two  orders  of  magnitude 
smaller.  The  large  value  of  r^,  with  its  corresponding 
value  of  x  =  0.463  at  one  centimeter,  has  placed  the  par¬ 
ticle  distribution  in  the  Ray leigh— Mie  transition  region. 
Note  also,  that  for  the  larger  value  of  refractive  index, 
the  ratios  of  absorption  to  scattering  coefficients  are 
lower  than  for  the  smaller  value  of  refractive  index. 

The  significance  of  these  results  is  that  cloud  emission 
can  be  accurately  calculated  in  a  straightforward  manner 
only  for  the  3  cm  case  with  the  postulated  particle  dis¬ 
tribution. 

Cloud  Albedo 

Under  the  assumptions  employed,  the  integrated  albedo 
for  single  scattering.  A,  for  the  model  cloud  is  qiven 
by 


so  that  at  T  +11  minutes 
o 


m  =  1 . 88-iO , 005  in  =  1. 63-iO.  005 


2.29x10' 


1 . 13xl0”2 


8.97x10' 


4,52x10* 
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and  at  T0  +  48  seconds 

m  «  1 . 86  -iO . 005  m  =  1  > 63-iO , 005 

Ak  =  2. 27x10“ 1  1. 12x10”* 

Ax  «  8. 93xl0”2  4.4x10“* 

The  integrated  albedo  provides  another  means  of  assessing 
both  the  emission  and  reflection  properties  of  the  cloud. 
For  example,  A  *  0.1  indicates  that  the  complete  radiative 
transfer  problem,  including  scattering,  absorption  and 
emission,  3hould  be  formulated  to  obtain  correct  results. 
The  amount  of  energy  from  ? lluminating  sources  scattered 
by  the  cloud  is  dependent  cn  source  strength  and  the  op¬ 
tical  thickness  and  albeac  of  the  cloud.  Thus,  while  the 
magnitude  of  the  albedo  alone  does  not  allow  determina¬ 
tion  of  the  likelihood  of  measurable  radar  return  from 
the  cloud,  the  backscattering  is  directly  proportional 
to  albedo.  Further  considerations  of  radar  backscatter¬ 
ing  are  described  in  Paragraph  7. 

d.  Cloud  Apparent  Temperatures 

It  has  been  concluded,  for  the  postulated  dust  particle 
distribution,  that  only  at  3  cm  can  it  be  assumed  that 
extinction  is,  for  practical  purposes,  entirely  by  ab¬ 
sorption.  Then,  at  this  wavelength,  if  it  is  further 
assumed  that  the  cloud  has  thermalized  at  temperature 
Tc,  the  cloud  apparent  temperature,  Tw,c,  would  be  given 
approximately  by  equation  (10).  Furthermore,  letting 
the  temperature  of  the  cloud  be  290°K,  the  cloud  apparent 
temperature  at  3  cm  for  cloud  parti  :les  with  the  refrac¬ 
tive  index  of  fine  sand  (m=l» 63-i0 . 005)  is 

Tbc  =  (1“.0987)  290  =  3.77°K 

The  total  apparent  temperature,  T^,  viewing  the  cloud 
against  the  sky  which  has  an  apparent  temperature,  T^, 
is  given  by 

Tb  “  Tbbt  +  ?bc  (13) 

while  the  change  in  apparent  temperature  from  the  sky 
background  to  occulted  background  is  given  by  equation 
(9).  For  an  elevation  angle  of  17°  (zenith  angle  =  73®), 
the  sky  brightness  would  be  approximately  14°K  (see  Figure 
7) .  Taking  this  value  as  the  average  sky  temperature  for 
calculation  purposes,  the  total  apparent  temperature  and 
change  in  apparent  temperature  for  the  postulated  dust 
cloud  T_  • r  45  seconds  axe: 
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Tg  =  14(0.987)  +  3.77  =  17.66K 

ATg  =  3.77  -  (0.013)14  =  3.59°K 

This  change  in  apparent  temperature  is  close  to  that 
actually  observed  during  Event  Dial  Pack,  ana  with  a 
slight  increase  in  mass  density  the  model  could  be  made 
to  match  the  value.  However,  without  comparing  the 
apparent  temperature  at  a  second  wavelength  and  estab¬ 
lishing  a  clpud  scattering  comparison  (e.g.,  degree  of 
polarization)  an  infinite  number  of  particle  distribu¬ 
tions  can  be  forced  to  satisfy  this  one  criteria. 

If  the  fact  that  the  albedo  at  1  cm  is  greater  than  C.l 
is  ignored,  and  the  apparent  temperature  and  change  in 
apparent  temperature  are  calculated  as  though  the  trans¬ 
mittance  were  a  measure  of  absorptive  loss,  the  corre¬ 
sponding  results  are 

=  49 . 5°K 

At£  =  9 . 5°  K 

The  ratio  of  the  changes  in  apparent  temperature  at  the 
two  wavelengths  is  2,64  as  compared  with  the  observed 
ratio  of  2.46.  While  this  ratio  is  important,  the  accu¬ 
racies  of  the  record  scalings  and  knowledge  of  which  30 
GHz  radiometer  channel  had  improper  gain  enter  into  the 
measured  ratio  (i.s.,  if  the  H  channel  was  assumed  to  be 
correct,  the  measured  ratio  would  be  1.5). 

The  magnitudes  of  the  preceding  results  are  representa¬ 
tive  of  those  actually  observed,  and  of  course,  the  ob¬ 
servations  influenced  certain  assumptions.  Besides  the 
assumption  that  the  Rayleigh  mechanism  holds  and  absorp¬ 
tion  uominates  over  scattering,  these  results  are  strong¬ 
ly  affected  by  the  choice  of  dielectric  properties  of  the 
scatterers,  their  size  distribution  and  their  number 
density. 

e.  Discussion 

The  length  and  detail  of  the  preceding  analysis  has  been 
intended  to  convey  the  thought  processes  as  well  as  the 
mechanics  of  arriving  at  a  solution.  This  approach  was 
taken  because  the  observational  data  suggested  that  ab¬ 
sorption  (insertion  loss)  was  the  do..  ’  nant  process,  i.e., 
the  measured  changes  in  apparent  temperature  caused  by 
the  uust  cloud  and  their  ratio  at  two  wavelengths  implied 
a  Rayleigh  mechanism.  Starting  with  cloud  particle  data 
taken  during  the  later  phase  of  cloud  evolution,  an 
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approximate  solution  for  an  early  stage  was  arrived  at 
in  two  steps.  This  solution  has  a  corresponding  particle 
size  distribution  which  has  a  significant  number  of  par¬ 
ticles  in  the  Rayleigh-Mie  transition  region  for  1  cm 
wavelength,  and  could  therefore  also  produce  a  degree  of 
polarization  at  this  wavelength.  It  must  be  emphasized 
though  that  this  solution  is  not  unique. 

Numerous  assumptions  were  made  throughout  the  analysis, 
many  of  them  due  to  the  lack  of  supporting  data,  but 
these  assumptions  had  no  drastic  consequence  because  of 
the  Rayleigh  behavior  of  the  real  cloud.  A  general  and 
complete  solution  which  includes  Mie  theory  and  formula¬ 
tion  of  the  radiative  transfer  problem  would  require 
supportive  data  in  lieu  of  certain  of  these  assumptions 
(such  as  homogeneity) . 

In  the  following  sections  the  results  of  Mie  calcula¬ 
tions  performed  as  a  check  on  the  preceding  Rayleigh 
approximations  and  to  extend  the  calculations  to  a  wider 
va  iety  of  particle  distributions  and  refractive  index, 
are  prcconted.  Tne  program  assembled  for  these  calcula¬ 
tions  will  form  the  nucleus  for  development  of  a  full 
general  solution. 


6.  MIE  THEORY  CALCULATIONS 

Even  though  the  particle  sizes  for  the  later  stage  cloud 
satisfy  the  Rayleigh  criterion,  those  in  the  earlier  phase  T0  + 

48  seconds  cloud  may  not.  Therefore,  full  Mie  theory  calculations 
have  been  made,  using  a  computer  program  for  which  the  accuracy 
has  been  verified  to  particles  of  size  parameter,  x,  greater  than 
5000,  Analysis  of  the  early  phase  cloud  necessitates  use  of  the 
program,  while  application  to  the  later  stage  cloud  provides  use¬ 
ful  checks  on  the  Rayleigh  approximations. 

a.  Description  of  Computer  Routines 

A  computer  program  was  written  in  FORTRAN  IV  to  calcu¬ 
late  the  scattering  and  absorption  properties  of  a  homo- 
ganeous  dust  cloud  made  up  of  spherical  particles  with 
an  arbitrary  size  distribution.  The  program  is  made  up 
of  several  separate  routines  which  will  be  briefly  de¬ 
scribed  here. 

A  routine  written  by  Dave,3  has  been  used  to  calculate 


3Dave,  J.  V.,  "Subroutines  for  Computing  Electromagnetic  Radia¬ 
tion  Scattering  by  a  Sphere,"  IBM  Scientific  Center,  Report  320- 
3237,  May,  1968. 
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the  Mie  scattering  parameters  for  each  particle  size  in 
the  dust  cloud  distribution.  Included  are  the  single 
particle  efficiency  factors  for  extinction  and  scatter¬ 
ing,  Qext(rfm)  and  Qsc(r»n')  from  which  the  absorption 
efficiency  factor  is  calculated  by  an  analogous  form  of 
equation  (7),  i.e., 

Qab  =  Q ext  “  Qsc  (14) 

In  addition,  four  elements  of  the  transformation  matrix, 
Van  de  Hulst,1*  are  determined  as  a  function  of  scatter¬ 
ing  angle,  and  from  these  the  intensity  and  degree  of 
polarization  of  scattered  radiation  are  calculated  for  a 
sphere  illuminated  by  an  unpolarized  beam  of  radiation. 

Once  the  scattering  and  extinction  efficiencies  have  been 
found  for  each  particle  size  in  the  cloud,  the  volume  co¬ 
efficients 


8[n(r),m]  =  tt  J  r  2n  (r )  Q  (r  .m)  dr  (15) 

r=o 


are  calculated  for  a  given  particle  density  distribution 
n(r).  In  the  program,  analytical  distributions  are 
numerically  integrated  using  Gauss-Legendre  quadrature, 
tabular  distributions  are  integrated  using  Lagrangian 
interpolation,  and  discrete  distributions  are  numerically 
summed. 

If  the  scattering  is  determined  to  be  negligible  (Bsc<< 
&bc)  with  respect  to  absorption,  and  assuming  a  homoge¬ 
neous  cloud,  the  apparent  cloud  temperature  is  calculated 
as  in  Paragraph  5.d,  That  is,  the  cloud  optical  thick¬ 
ness  is  calculated  by  equation  (4) ,  cloud  transmit ta nee 
by  equation  (5),  cloud  emission  by  equation  (10)  and  the 
apparent  temperature  of  the  cloud  by  equation  (13). 

■  The  program  also  provides  other  measures  of  the  cloud, 
such  as  mass  of  particles  per  unit  volume,  equivalent 
volume  radius  and  equivalent  volume-squared  radius  for 
the  given  particle  distribution. 

b.  Late  Phase  (T0  +  11  Minues)  Cloud  Results 

Using  the  discrete  distribution  for  pass  #1,  the  scatter- 


4 

Van  de  Hul3t,  H.  C. ,  "Light  Scattering  by  Small  Particles," 
John  Wiley  and  Sons,  Inc.,  1957. 
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ing  and  brightness  parameters  were  calculated  for  a  di¬ 
electric  constant  ot  2.3-10.023  and  for  the  coarse  and 
fine  sand  dielectric  constants,  with  imaginary  parts  of 
0.001,  0.005,  and  0.01.  Table  IV  summarizes  the  results 
of  that  analysis. 


c.  Early  Phase  Cloud  Model  and  Results 


With  the  exception  of  the  observed  brightness  temperature 
changes  at  T0  +  48  seconds,  there  is  very  little  data 
from  which  a  particle  distribution  at  the  early  stage 
can  be  derived.  Because  of  this  lack  of  experimental 
information,  another  approach  has  been  taken  to  formu¬ 
late  an  early  stage  cloud  model. 


An  analytical  distribution  often  used  to  represent  par¬ 
ticulate  clouds  is  the  generalized  gamma  function5 

n(r)  =  araexp{-br*0  (16) 


with 


a  = 


,N 


b  (a+1 )  /y 


v  ia+l> 


=  3  r-Y 
Y  C 


where  rc  is  the  mode  radius,  and  N  is  a  density  para¬ 
meter.  By  varying  the  constants  a  and  y,  this  function 
can  be  made  to  fit  many  observed  cloud  particle  distri¬ 
butions,  over  a  large  range  of  particle  sizes.  For  the 
values  of  r  <  rc,  the  shape  of  the  distribution  is  deter¬ 
mined  primarily  by  the  value  of  a,  while  that  for  r >  rc 
is  determined  largely  by  y. 


Using  the  gamma  function  as  an  envelope  for  the  discrete 
distribution  of  pass  #1,  the  values  of  a  =  1.2  and  -y  = 
1.2  have  been  determined  (Figure  12). 

In  formulating  the  early  phase  cloud,  there  is  good 
reason  to  believe  that  only  the  density,  and  not  the 
shape  of  the  distribution  for  smaller  particles  will 
change  from  T0  +  48  seconds  to  T0  +  11  minutes.  On  the 
other  hand,  the  shape  of  the  distribution  for  larger 
particles  will  not  yet  be  stable  at  T0  +  48  seconds,  and 
will  change  as  the  larger  particles  settle  out. 


sDeirmendjian,  D. ,  "Complete  Microwave  Scattering  and  Extinction 
Properties  of  Pclydisperssd  Cloud  and  Pair.  Elements. "  The  Rand 
Corporation,  Report  R-422-PR,  December  1963. 
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Figure  12  -  Late  Phase  Cloud  Distribution 
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With  this  in  mind,  the  early  phase  cloud  distribution 
was  determined  using  an  iterative  process:  a  was  held 
constant  at  1.2,  N  and  rc  were  varied  over  a  range  of 
values  consistent  with  the  extrapolated  mass  and  si2e 
data,  (Figure  10),  and  y  was  varied  between  1  and  2. 
Figure  13  shows  the  distributions  for  four  dielectric 
constants  which  give  the  changes  in  brightness  temper¬ 
ature,  Tfc,  from  background  to  occulted,  as  shown  in 
Table  V.  The  corresponding  volume  coefficients  are  also 
shown  in  Table  V. 

Comparing  the  absorption  and  scattering  coefficients  at 
1  cm,  the  albedo  is  greater  than  0.1  in  three  cases,  im¬ 
plying  that  the  full  radiative  transfer  equation  should 
be  solved  to  determine  a  uore  accurate  brightness  tem¬ 
perature.  At  2.94  cm,  the  scattering  coefficient  is 
smaller  than  the  extinction  coefficient  by  at  least  an 
order  of  magnitude,  and  the  self emission  temperature  of 
the  cloud  as  calculated  is  valid. 

d.  Comparison  of  Mie  and  Rayleigh  Results 

Excellent  agreement  can  be  seen  between  the  volume 
scattering  and  extinction  coefficients  calculated  for 
the  late  stage  cloud  model  using  the  full  Mie  solution 
(Table  IV) ,  and  those  calculated  using  the  Rayleigh  ap¬ 
proximations.  In  all  cases,  the  approximation  error  is 
less  than  a  few  percent,  and  the  assumption  of  an  equiv¬ 
alent  monodispersed  medium  with  particles  of  an  effective 
volume  radius  and  volume-squared  radius  is  certainly 
valid. 

On  the  other  hand,  the  early  phase  cloud  particle  size 
distribution  model  which  has  been  used  shows  evidence 
that  the  scattering  process  is  becoming  non-Rayleigh.  A 
measure  of  this  departure  from  Rayleigh  scattering  is  the 
ratio  of  extinction  coefficients  *  For  the  late 

phase  cloud  distribution 

8ex/&ex  2.98 

The  reason  the  ratio  is  somewhat  greater  than  the  ex¬ 
pected  value  of  2,94  is  probably  due  to  the  presence  of 
the  few  large  particles  which  approach  the  Rayleigh  size 
limit.  For  the  early  phase  cloud  motel  distributions 
(Figure  13) ,  the  ratio  has  significantly  increased  to 

Sex/Sex  3.5 

This  is  ar.  expected  result  since,  if  tne  Mie  process  is 
pr cuuiuuia te  at  Luc  hiy/ttsi  frequency,  the  extinction 
will  be  significantly  larger  then  that  predicted  by  the 
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Particle  Density  (cm" 


Particle  Radius  (u) 


Figure  13  -  Early  Phase  Cloud  Distributions 
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Rayleigh  approximations . 

Similar  analysis  using  the  absorption  coefficients  for 
the  later  phase  cloud  shows  that 


e 


k 

ab 


2.95 


whereas  the  early  stage  distributions  all  yield 

e£b  '  3-<>5 

This  is  further  evidence  of  the  non-Rayleigh  nature  of 
the  cloud. 


Viewing  these  results,  it  is  now  apparent  that  the  ratio 
of  2.64  found  in  Paragraph  5.d  for  the  ratio  of  the 
changes  in  apparent  temperature  at  1  cm  and  3  cm  is  low 
for  the  reason  that  the  Rayleigh  approximations  are  in 
error  for  particles  (re,  rw)  which  violates  the  Ray¬ 
leigh  criteria  x  <  0.1.  At  the  same  time,  the  question 
as  to  why  the  ratio  obtained  from  the  Event  Dial  Pack 
measurements  was  even  lower  is  raised.  The  answer  is 
not  clear,  but  the  alternatives  are  calibration  errors 
or  properties  of  the  dust  cloud  not  accounted  for  (e.g., 
scattering  or  structure) . 


e.  Columnar  Dust  Content 


The  quantity  of  layered  sand  required  to  cause  comparable 
increases  in  the  30  GHz  apparent  temperatures  of  the 
zenith  sky  to  those  observed  on  the  Event  Dial  Pack  cloud 
are  0.387  gm  cm-2  and  1.06  gm  cm”2  for  coarse  and  fine 
sand,  respectively,  (see  Figure  8) .  The  approximate 
values  obtained  from  the  computer  models  are  found  by 
dividing  the  mass  of  particles  per  unit  volume  by  the 
particle  equivalent  volume  per  unit  volume.  The  results 
for  coarse  and  fine  sand  are  25  qm  cm"  and  27  gm  cm"  , 
respectively.  These  results  are  tractable  because  bulk 
material  attenuates  more  than  the  same  material  in  sus¬ 
pension.  Furthermore,  the  ratio  of  bulk  extinction  to 
particulate  extinction  in  the  Rayleigh  region  can  be 
several  orders  of  magnitude.  This  ratio  for  the  same 
material  at  wavelengths  that  place  the  particulates  in 
the  Mie  region  can  approach  unity.  Therefore,  for  par¬ 
ticle  size  distributions  in  the  Rayleigh-Mie  transition 
region,  an  order  of  magnitude  seems  reasonable. 


7 .  RADAR  BACKSCATTER 

Tilts  i  dual  ie  LuiTi  fiJOIu  tTiS  uU5t 

using : 


—  1  —  .  .  J  ---  U  1 
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P  A28<J>£ 

75lV-  escexP(-2eexL>  (17) 

received  (backscattered)  and  transmitted 
powers 

antenna  apartural  area 

horizontal  and  vertical  antenna  bandwidths 
radar  pulse  length 

range  to  the  centroid  of  the  scattering 
volume 

volume  scattering  coefficient  of  the  cloud 

range  to  the  centroid  of  the  scattering  vol¬ 
ume  from  the  leading  edge  of  the  cloud 

The  following  characteristics  of  the  LAPQ-1  radar,  used  in  the 
Event  Dial  Pack,  are  known: 

Peak  power  =  IOC  kw 

Wavelength  ■=  .858  cm 

Pulse  length  -•  1  sec 

Minimum  detectable  signal  =-101  db 

Antenna  gain  =  47  db 

Antenna  beamwidth  =  0.75  degrees 

Applying  equation  (17)  for  a  radar  range  of  3000  m,  the  approxi¬ 
mate  radar  return  from  the  dust  cloud  is  given  by 

pr  =  .92  6  (18) 
r  sc 

Evaluating  this  equation  for  a  cloud  with  particle  refractive 
index: 


where: 


4  r,t 

Ag 

6,4 

£ 

R 


e 


sc 


m  -  1.63  -  iO.Ol 

and  the  volume  scattering  coefficients  corresponding  to  the  late 
phase  cloud  data,  and  the  postulated  early  phase  cloud  model 
(Paragraph  6.c); 
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the 


3 

sc 


1. 53x10  1 3  cm  1 t  T0  + 
8.44xl0”e  cm’1)  Tq  + 


11  minutes 
48  seconds 


radar  returns  are: 


Pr  =  1.41xlo”13  watts ;  Tq  +  11  minutes 

which  is  approximately  -98  dbm  during  the  dust  cloud  dissipation 
phase,  and 


P_  =  7.76x10  8  watts;  T  +  4 8  seconds 

or  approximately  -41  dkm  during  the  growth  phase.  This  indicates 
that  for  an  LAPQ-1  type  radar,  the  cloud  should  be  readily  detect¬ 
able  past  Tc  +  11  minutes.  One  factor  not  taken  into  account  is 
a  possible  change  of  radar  range  as  the  cloud  drifts  from  ground 
zero.  A  range  increase  of  10  times  (i.e.,  R  =  30,000  m)  would 
have  made  the  cloud  undetectable  at  T0  +  11  minutes. 

For  a  radar  operating  at  3  cm  with  the  LAPQ-1  system  char¬ 
acteristics  ,  the  radar  return  is  given  by 


Pr  =  10.1  83C  (19) 

With  a  particle  refractive  index  m  a  1.63-i.Ol,  the  scattering 
coefficients  are: 


I"  1.11x10  15  cm  1 ;  T 0  +  3.1  minutes 

6  3q  I  — in  — l 

j  6.07x10  cm  ;  Tc  +  48  seconds 

The  radar  returns  are  than 

Pr  =  1.12x10  14  watts;  TQ  +  11  minutes 

or  about  -103  -bm  for  the  late  stage  cloud,  and 

Pr  «  6.12x10  5  watts;  T0  +  48  seconds 

or  about  -52  dbm.  The  early  phase  cloud  is  still  readily  detect¬ 
able,  while  the  late  stage  cloud  would  probably  not  be  seen. 


8.  POLARIZATION  CONSIDERATIONS 

The  approximate  methods  used  in  the  analysis  up  to  this 
point  tacitly  assume  a  spherical  homogeneous  cloud  composed  of 
spherical  particles.  If  it  is  additionally  assumed  that  single 
scattering  governed  by  the  Rayleigh  laws  applies,  and  that  the 
cloud  is  optically  thin,  polarized  emission  would  not  be  ob¬ 
served.  Even  if  the  particles  were  non-spherical,  their  random 
orientation  in  the  cloud  should  eliminate  any  polarized  emission. 
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Since  the  cloud  was  non-homogeneous  and  not  spherical,  it  is 
possible  to  imagine  a  situation  in  which,  although  the  cloud  is 
optically  thin,  the  radiometer  could  be  viewing  stratification 
within  the  cloud  at  near  gracing  angles.  Depending  on  the  struc¬ 
ture  and  the  density  gradients  which  the  radiometer  is  viewing, 
some  polarized  emission  could  be  observed.  However,  such  an 
occurance  would  be  evidenced  by  an  irregular  brightness  curve,  as 
the  cloud  passed  through  the  radiometer  beam,  while  the  Event  Dial 
Pack  record  is  relatively  smooth.  This  leaves  the  most  likely 
source  of  polarized  energy,  if  observed,  the  scattering  of  radia¬ 
tion  which  illuminates  the  cloud. 

The  sources  of  cloud  illumination  to  be  considered  are  the 
sun,  the  atmosphere,  and  the  ground,  since  the  radars  used  to 
probe  the  cloud  have  been  treated  separately  in  Paragraph  7. 
Polarized  or  unpolarized  illuminating  radiation  can  give  rise  to 
measured  polirization  if  the  source  strengths  are  sufficient  for 
the  scattered  energy  to  be  detectable.  Highly  polarized  incident 
radiation  will  in  general  be  somewhat  depolarized,  while  unpolar¬ 
ized  incident  energy  will  exhibit  some  degree  of  polarization 
when  scattered  from  a  cloud  such  as  is  being  discussed  here. 

The  direct  sun  is  the  easiest  case  to  treat,  and  its  treat¬ 
ment  will  give  some  insight  to  the  problem.  At  1  cm  wavelength 
the  sun  can  be  thought  of  as  a  distant  point  source  with  a 
brightness  temperature  of  about  7000°K.  Referring  to  Figure  6, 
which  shows  the  experiment  gecmetry,  it  is  seen  that  the  half 
power  radiometer  beam  defines  an  approximately  cylindrical  volume 
within  the  cloud.  The  particles  contained  in  this  volume  are 
assumed  to  be  the  only  active  scatterers.  The  dimensions  of  the 
cylinder  are  1000m  length  (the  width  of  the  cloud  at  TQ  +  48  sec¬ 
onds)  by  113m  mean  diameter. 

To  determine  if  scattered  solar  radiation  could  contribute 
to  the  measured  brightness  of  the  dust  cloud,  a  best  possible 
case  is  considered.  Let  the  scattering  angle  be  90  degrees  so 
that  all  elements  of  the  scattering  phase  matrix  are  zero  except 
Pj..  The  power  density  of  the  solar  radiation,  C,  incident  on  the 
dust  cloud  at  1  cm  wavelength  is  of  the  order  of  10“25  watts 
cm"2Hz“l.  Under  the  single  scattering  assumption,  each  cubic 
centimeter  of  the  cylindrical  volume  will  be  illuminated  by  radia¬ 
tion  of  this  intensity,  and  each  will  scatter  an  amount  l/2BscPiS 
toward  the  radiometer  (ignoring  extinction  of  the  incident  or 
scattered  power  by  the  cloud).  The  power  density  of  the  scattered 
solar  flux,  Ss ,  at  the  radiometer  is 

8  P  SV 

Ss  =  1/2  ~~~ —  (20) 

4ttR2 

where  V  is  the  volume  of  the  scattering  cylinder  and  R  is  the 
range  of  the  radiometer  from  the  centroid  of  the  cloud.  Evalua- 
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tinq  equation  (20)  at  T0  +  48  seconds  and  X  =  1  cm  for  the  postu¬ 
lated  early  phase  cloud  particle  distribution  (see  Paragraph  6.b) 
with  refractive  index  m  =  1.63-i0.01  and  6SC  =  4.56x10“®  cm”1  (see 
Paragraph  6.b),  gives 


3.74x10” 


watts  cm' 


The  power  density  at  the  radiometer  corresponding  to  the  apparent 
temperature  of  the  cloud  is  ^iven  by 

2irkTb  2 

Sc  .  -p—  (5>  (21) 

where  r  is  the  radius  of  the  scattering  cylinder  defined  by  the 
antenna  beamwidth,  and  k  is  Boltzmann's  constant.  Evaluation  of 
this  equation  gives 

S  =  l.SixiO-24  watts  cm”2  Hz"1 

c 

which  is  eight  orders  of  magnitude  greater  than  the  scattered  solar 
flux.  Therefore,  it  can  be  concluded  that  scattered  solar  energy 
can  be  of  no  consequence. 

Consideration  of  radiation  from  the  ground  is  somewhat  more 
complicated.  To  make  the  analysis  tractable,  the  ground  is 
treated  as  an  idealized  plane  Lambertian  surface  of  uniform  tem¬ 
perature,  Te.  The  power  density  of  the  radiation  from  the  ground 
at  the  centroid  of  the  cloud  is  approximately 

c  =  2rkT 


where  it  has  been  assumed  that  the  effective  surface  is  contained 
in  a  solid  angel  of  tt  steradians  as  a  consequence  of  the  Lamber¬ 
tian  properties  of  the  ground.  Evaluation  of  the  equation  (22) 
for  a  surface  temperature  of  300°K  gives 

S£  =  2.60xl0“20  watts  cm"2  Hz"1 

which  yields  the  corresponding  power  density  at  the  radiometer; 

s|  <  6.48xl0'"27  watts  cm”2  Hz-1 

Note  that  this  value  has  been  explicitly  stated  as  "less  than" 
because  of  numerous  optiirium  assumptions.  Comparing  this  value 
with  the  power  density  corresponding  to  the  cloud  apparent  tem¬ 
perature,  it  is  seen  to  be  more  than  two  orders  of  magnitude 
smaller.  This  leaves  the  terrestrial  radiation  scattered  by  the 
cloud  at  the  detection  threshold  of  the  radiometer,  with  negli¬ 
gible  contribution  to  the  total  signal. 

Atmospheric  radiation  presents  a  similarly  difficult  radia- 
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tive  problem,  in  which  the  cloud  is  immersed  in  the  source.  How~ 
ever,  the  radiation  temperature  of  the  atmosphere  averaged  over 
4^  steradians  is  less  than  100°K  at  1  cm,  which  produces  a 
scattered  power  density  at  the  radiometer  of  about  the  same  level 
as  the  terrestrial  radiation.  Therefore,  it  can  be  concluded  that 
radiation  from  the  atmosphere  can  also  be  neglected. 

The  treatment  of  these  various  sources  was  based  on  the  selec 
tion  of  one  particle  size  distribution  that  led  to  reasonable 
agreement  with  observations.  While  this  does  not  allow  unique 
assessment  of  the  cloud  emission  and  scattering,  the  results  seem 
to  be  compatible  with  the  limited  data  available. 


SECTION  III 


SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


In  the  preceding  analyses,  the  problem  of  microwave  emission 
and  scattering  by  a  tenuous  cloud  of  earth  debris  (called  dust 
throughout  the  report)  has  been  considered  in  an  approximate  way. 
First,  using  debris  sample  data  acquired  with  an  airborne  par¬ 
ticle  sampler  and  particle  refractive  index  values  the  volume  ex¬ 
tinction,  scattering  and  absorption  coefficients  were  calculated 
by  Rayleigh  formulas.  In  addition,  optical  thickness,  trans¬ 
mittance  and  cloud  albedo  were  calculated.  Second,  an  extrapola¬ 
tion  to  an  early  time  (T0  +  48  seconds)  was  made  by  increasing 
the  mass  density  by  a  factor  of  4.4x10**  times.  Again,  the  radia¬ 
tion  transfer  parameters  were  calculated  by  the  Rayleigh  approxi¬ 
mations.  This  extrapolation  led  to  cloud  transmittances  and 
cloud  brightness  temperatures  close  to  those  actually  observed  by 
the  10.2  and  30  GHz  radiometers  during  Event  Dial  Pack. 

As  a  check  on  the  Rayleigh  results  and  to  extend  the  calcu¬ 
lations  to  a  wider  range  of  particle  size  distributions  and  re¬ 
fractive  indices,  a  computer  program  was  written  utilizing  gen¬ 
eralized  Mie  scattering  routines.  With  a  version  of  this  program, 
the  validity  of  the  Rayleign  approximations  for  the  late  stage 
(T0  +  11  minutes)  cloud  was  determined.  In  addition,  a  dust 
cloud  particle  distribution  model  was  postulated,  and  using  this, 
a  combination  of  mass  density,  particle  density  and  refractive 
index  were  iteratively  selected  to  yield  results  close  tc  the 
actual  observations  of  the  early  stage  (Tc  ->•  4  8  seconds)  cloud. 

Considerations  of  radar  detectability  using  computed  values 
of  Mie  scattering  coefficients  suggest  rhat  the  c  loud  should 
have  been  readily  detectable  in  its  early  stages  by  both  3  and 
0.86  cm  radars,  and  therefore  information  about  its  structure 
could  be  determined  by  such  means.  In  the  later  phase  of  cloud 
evolution,  T  >  T0  +  11  minutes,  the  cloud  should  still  have  been 
giving  a  usable  return  of  0.86  cm  buc  no'-,  at  3  cm.  It  is  not 
known  whether  this  is  the  case,  but  the  unique  determination  of 
early  phase  cloud  properties  hi  ages  on  this  type  of  data. 

Using  the  same  model  as  was  used  to  assess  radar  return,  the 
scattering  of  radiation  from  the  sun,  ground  and  atmosphere  toward 
the  radiometers  was  treated.  None  of  these  sources  was  found  to 
contribute  measurable  amounts  of  radiation,  and  therefore  do  not 
contribute  any  significant  polarized  energ  '  for  the  postulated 
early  stage  cloud  model. 

The  conclusions  of  this  report  are  as  follows: 

1.  The  late  phase  cloud.  T  >  Tc  +  11  minutes,  would  have 
been  transparent  to  th  10,2  and  30  GHz  radiometers. 
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2.  There  are  many  models  which  approximately  predict  the  ob¬ 
served  apparent  temperatures  of  the  early  phase  cloud 

(T0  +  48  seconds )  . 

3.  The  controlled  measurements  data  show  that  polarization 
was  insignificant.  The  derived  models  are  consistent 
with  this  conclusion. 

4.  Approximately  25  gm  cm-2  of  debris  was  required  to  cause 
the  observed  increases  in  apparent  temperature. 

5.  Strong  radar  returns  should  have  been  observed  at  0.86 
and  3  cm  from  the  early  phase  cloud,  but  the  late  phase 
cloud  should  have  been  marginally  detectable  only  at 
0.96  cm. 

6.  Insertion  loss  and  transmission  loss  cannot  be  absolutely 
calculated  because  of  cloud  backgrounds  and  change  in 
cloud  thermal  temperature,  cloud  transmittance  is  a 
better  descriptive  property. 

If  reliable  polarization  and  radar  backscatter  data  had  been 
available,  it  would  have  been  possible  to  derive  a  unique  model  of 
a  homogeneous  spherical  debris  cloud.  To  handle  a  structured 
cloud  or  one  in  which  scattering  was  an  important  process  would 
have  requirt d  formulation  of  the  radiative  transfer  problem,  a 
task  oeyond  the  scope  of  this  study.  As  a  spherical,  homogeneous 
cloud  was  assumed,  the  complete  scattering  by  the  cloud  would  be 
given  classical  scattering  patterns. 

The  value  of  the  radiometers  in  experiments  such  - s  Event 
Dial  Pack  should  not  be  underestimated,  for  they  provide  informa¬ 
tion  complementary  to  that  obtained  by  radar  techniques.  Inter¬ 
pretation  of  these  two  types  of  data  jointly  should  provide 
details  on  the  structure  of  such  a  cloud,  as  well  as  properties 
of  the  particles  of  which  it  is  composed.  Whereas  a  radar  can 
supply  scattering  and  extinction  data,  a  radiometer  can  supply 
absorption  and  extinction  data,  which  together  completely  define 
the  radiative  properties  of  the  cloud.  It  is  therefore  recom¬ 
mended  for  future  experiments  of  this  type,  that 

1.  Three  radiometers,  one  each  at  wavelengths  of  3,  1.25 
and  1  cm,  be  installed  at  the  test  location.  The  1.25 
radiometer  would  be  used  to  provide  data  on  water  vapor 
and  cloud  thermalization  as  well  as  providing  a  addi¬ 
tional  constraint  for  iterative  solution. 

2.  These  radiometers  should  be  absolutely  calibrated  to 
give  the  total  apparent  temperature  as  well  as  the 
change  a n  apparent  temperature 

3.  These  radiometers  should  have  the  capability  of  scanning 
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whenever  the  cloud  has  dimensions  large  compared  with 
the  radiometer  beamwidth. 

4.  These  radiometers  should  track  the  cloud  as  it  ascends 
and  then  drifts  from  ground  zero. 

5.  Synchronizing  information  between  radar  positioning  and 
radiometer  positioning  be  recorded  for  analysis  purposes. 

In  addition  to  the  refined  radiometric  systems,  a  complete 
general  computer  program  which  includes  the  radiation  transfer 
formulation  for  an  arbitrarily  shaped  cloud  should  be  developed. 
Some  of  the  basic  elements  for  such  a  computer  model  are  already 
available,  including  the  scattering  routines  utilized  in  this 
study  and  radiation  transfer  models  utilized  in  past  studies. 

With  a  complete  model  available,  coordination  among  experimenters 
and  proper  data  recording  and  handling  could  assure  rapid  accom¬ 
plishment  of  the  dust  cloud  analysis. 


